Pattern formation inside bacteria: fluctuations due to low copy number of proteins 
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We examine fluctuation effects due to the low copy number of proteins involved in pattern-forming 
dynamics within a bacterium. We focus on a stochastic model of the oscillating MinCDE protein 
system regulating accurate cell division in E. coli. We find that, for some parameter regions, the 
protein concentrations are low enough that fluctuations are essential for the generation of patterns. 
We also examine the role of fluctuations in constraining protein concentration levels. 
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In recent years, dramatic experimental progress has 
been made in resolving the subcellular localization of bac- 
terial proteins. Often the proteins form self-organized, 
spatially inhomogeneous patterns that can involve coher- 
ent spatiotemporal oscillations. These self-organized pat- 
terns are vital for accurate cell division i^^^^,^ Jj. 
Understanding these patterns promises to reveal new 
mechanisms for the generation of subcellular bacterial 
structure. In this letter, we examine, for the first 
time, the impact of fluctuations on these patterns, fo- 
cusing on the oscillatory MinCDE system in E. coli 

QiiiiiQ. 

Each E. coli cell divides roughly every hour, first repli- 
cating its DNA into two nucleoids, and then dividing at 
midcell into two daughter cells. If division is not targeted 
accurately to midcell then DNA will not be distributed 
to both daughter cells, resulting in unviable anucleate 
"minicells" . Division is initiated by a polymeric ring of 
the protein FtsZ, which forms on the inner side of the 
cytoplasmic membrane. Precise positioning of the FtsZ 
ring to midcell is controlled both by the inhibiting ef- 
fect of the nucleoids ("nucleoid occlusion"), and by the 
MinCDE system of proteins "s^ . MinC inhibits the forma- 
tion of the FtsZ ring, and is recruited to the membrane 
by MinD. MinE is also recruited to the membrane by 
MinD, where it forms a dynamical ring structure that 
drives pole to pole oscillations of MinC and MinD, with 
a period of about a minute. The oscillations lead to 
a time-averaged midcell MinC concentration minimum 
0. As MinC inhibits FtsZ ring formation, the FtsZ can 
only assemble near the cell midplane. The self-organized 
protein patterns are thus used by the cell to obtain posi- 
tional information without stationary positional markers 

mniiii 

Very recently, the physical principles behind the 
MinCDE protein patterns have been explored through 
reaction-diffusion equations 0, 0, 0|. The resulting 
equations represent protein diffusion, both in the cyto- 
plasm and along the membrane, and also protein bind- 
ing/unbinding from the cytoplasmic membrane. The 
slow membrane diffusion used in these models, relative to 
the cytoplasmic diffusion, results in Turing-like (Hopf) 



instabilities that spontaneously generate oscillatory pat- 
terns, in good agreement with experiment |9l llOl lll|. 

Bacterial proteins are, however, typically present in 
low numbers within the cell. This induces large fluctu- 
ations [l^ IT^ which have not been considered in pre- 
vious pattern-forming models |^ Hol . In E. coli, a 
recent assay put the copy numbers for MinD and MinE 
at 2000 and 1400 respectively We have therefore 

investigated the role of fluctuations in a discrete particle 
model of the E. coli MinCDE system where each protein 
molecule is explicitly tracked. This allows for a full anal- 
ysis of the fluctuations of both reactions and diffusion. 
Although the effects of noise have been studied in subcel- 
lular models without spatial dynami cs llSllig. and also 
in some spatially extended patterns |l7l llSl llol I20I I2H , 
the effects of fluctuations on subcellular positional infor- 
mation, and the constraints on protein concentrations 
due to fluctuations, have not previously been studied. 
We find that: 

(i) Even at surprisingly low concentrations, the noise 
does not destroy the oscillatory dynamics and indeed can 
be vital for generating patterns in regions of parame- 
ter space where the equivalent deterministic dynamics 
decays away. A bacterium can thus exploit low copy 
number fluctuations to produce stable, self-organized 
patterns. This result likely applies to any stochastic 
reaction-diffusion model of pattern formation involving 
sufficiently few protein copies (i.e. thousands or fewer). 

(ii) We find evidence that the cell employs sufficient copy 
numbers of the MinCDE proteins to ensure reliable mid- 
cell division; using substantially fewer copies degrades ac- 
curate positioning of the midcell MinCD minimum, using 
more does not lead to significantly improved accuracy. 

Stochastic model. We begin by introducing the 
stochastic model for the MinCDE dynamics, based on 
our deterministic model of Ref. We employ a Id 
discrete particle model, where the particles hop between 
lattice sites and where the full fluctuation effects are in- 
trinsically included by discrete particles. As in earlier 
models [3, E, 0| , the MinC dynamics is omitted, since 
it is known to closely follow the MinD dynamics. The 
occupancy at site i is JT-j'^, with j — {D,d,E,e} rep- 
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resenting cytoplasmic MinD, membrane MinD, cytoplas- 
mic MinE, and membrane MinE, respectively. Each pro- 
tein molecule is represented as a particle which, at each 
timestep, may hop with equal probability _Dj Ai/(Aa;)-^, 
where Ai is the time increment, to one of its neighboring 
sites at a; ^ a; ± Aa; (except for the boundary sites at 
either end, where hard wall boundary conditions are im- 
posed). At site i, the following reactions may occur, the 
first being for each D particle, and then for each d, E, e 
particle respectively: 
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These reactions are the stochastic analogs of the reac- 
tion processes used in our deterministic partial differen- 
tial equation model 0. The di term describes sponta- 
neous membrane association of MinD, whereas the (72 
term describes ejection of MinD from the membrane by 
membrane-bound MinE. Similarly, the (74 term describes 
spontaneous membrane disassociation of MinE, whereas 
the tT3 term describes recruitment of MinE to the mem- 
brane by cytoplasmic MinD. The a'l, a'^ "suppression" 
terms correspond to membrane MinE suppressing the 
binding of MinD to the membrane, and to cytoplasmic 
MinD suppressing the unbinding of membrane MinE. 

ATP dynamics. We next address the question of 
how ATP dynamics fits into the model. Experimentally, 
MinD is an ATPase HI and it is the MinD- ATP com- 
plex that binds to the membrane, whereas the release 
of MinD back into the cytoplasm requires MinE-induced 
ATP hydrolysis Our model assumes that, following 
ATP hydrolysis and release of MinD into the cytoplasm, 
nucleotide exchange is sufhciently rapid to allow for mem- 
brane reattachment of MinD-ATP almost immediately. 
As a result, we only model MinD-ATP in the cytoplasm 
(n^^) and on the membrane (n^^). Like "actin tread- 
milling" in eukaryotic cells, the ATP driven binding and 
unbinding of MinD allows for a cyclic MinCDE pattern 
to be maintained with only low levels of protein synthe- 
sis: this makes the pattern-forming dynamics extremely 
energy efficient. This recycling is also consistent with 
experiments where protein synthesis was blocked, but 
where the MinCDE oscillations were observed to continue 
unaffected 0. Consequently, the above model (similar 
to that of Ref. |0) does not include protein synthesis or 
degradation. These processes occur, but on longer time 
scales than the relatively rapid MinCDE oscillations. 

Simulations. In our stochastic model simulations, we 
use time and spatial increments Ax = 0.02 fj,m, At = 2 x 
10~^ s, so that 100 lattice sites model a 2 /im bacterium. 
We use: Do = 0.28 ^J,m^ s-\ Dd = 0.003 fim^ s-\ 



De = 0.6 /iTO^ s^i, De = 0.006 Aim2 s'^, ai = 20 s'^, 
and (T4 = 0.8 s~^. Note that the membrane diffusion 
constants are much smaller than those in the cytoplasm; 
this agrees with recent data indicating that MinD may 
polymerize on the membrane |23L l24j . For the remaining 
variables of the model, we focus on 4 representative pa- 
rameter sets shown in Table I, where we define N as the 
total number of MinD proteins, equal to the total num- 
ber of MinE proteins [2^. However, we emphasize that 
our results for the oscillatory behavior observed below 
are typical for large regions of parameter space. Initially, 
MinD and MinE particles are randomly distributed on 
the membrane and in the cytoplasm. Equal numbers of 
proteins are used since "wild-type" oscillations are ob- 
served when both proteins are equally expressed on plas- 
mids (this is consistent, within experimental uncertain- 
ties, with the earlier quoted MinDE assay 

Fluctuation driven instability. Using the parameters 
in Table I, we find that the presence of noise is vital for 
the oscillations to persist. This is shown in Fig. 1, where 
the ratio of the average MinD density in the right-hand 
30% of the cell to that in the left-hand 30% of the cell 
is plotted as a function of time for the stochastic model 
(at TV = 200 and N = 1500) and for the determinis- 
tic model with equivalent parameters (i.e. with the 
above reaction probabilities directly transformed into de- 
terministic reaction rates). In both cases, for the deter- 
ministic model, the protein concentrations rapidly decay 
away to the homogeneous steady state (in agreement with 
linear stability analysis Q), whereas regular oscillations 
continue for the stochastic model. Hence the average 
behavior of the stochastic model is clearly not describ- 
able using the naively equivalent deterministic model. 
To investigate this issue in more detail, we have exam- 
ined how steady-state/oscillation bifurcations in the de- 
terministic model are altered in the stochastic model. 
As a representative example, at iV = 1500, using the 
above parameter set in the stochastic model, but with 
Dd = 0.001 fim'^ s-\ De = 0.003 ^im^ 5-1, a[ ^ 0.2 and 
varying (T4, we find that the transition from oscillatory 
to steady-state behavior is reduced by around 40% from 
(T4 — 0.63 (deterministic) to (T4 « 0.39 s^^ (stochas- 
tic), as determined using the end to end MinD ratio. 
The noise does smear out the transition somewhat in the 
stochastic model (the transition at (T4 w 0.39 s^^ has 
width ±0.05 s~^), but this effect is rather small. Hence 
this smearing out cannot account for the large regions of 
parameter space where oscillations occur in the stochas- 
tic, but not the deterministic, model. Rather, we have a 
fluctuation driven instability, where the noise has shifted 
the location of the transition, thereby promoting oscilla- 
tions in large regions of parameter space where it would 
be forbidden in the equivalent deterministic model (see 
Fig. 1). Cells can in principle exploit low copy num- 
ber fluctuations to generate pattern-forming dynamics. 
The oscillations continue down to very low concentrations 
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TABLE I: Reaction rate parameter values. 

{N = 200), underlining the robustness of the dynamics. 

The above examples show that the fluctuations are 
often essential for pattern formation. However, it also 
possible for the stochastic and deterministic models with 
equivalent parameters both to generate oscillations. 
Hence we cannot definitively conclude that fluctuations 
are essential for the MinCDE oscillations. Nevertheless, 
our analysis does show conclusively that cells can exploit 
low copy number fluctuations for the generation of dy- 
namical subcellular structure. 

Effect of fluctuations on the midcell MinD minimum. 
In Fig. 2, we plot the MinD and MinE concentration 
profiles for N — 200 and N = 1500, showing their av- 
erages over 160 and 110 successive cycles, respectively, 
and also for 4 individual data sets each, where each set 
is averaged over only an individual oscillation cycle. For 
the long time data sets, we find that the midcell MinD 
concentration minimum (and a MinE concentration max- 
imum) are still robustly reproduced even in the presence 
of noise. However, as can be seen from the data aver- 
aged over individual cycles, the fluctuations around this 
average can be very large for small N. In Fig. 3 we 
show histograms of the position of the MinD concentra- 
tion minimum, where each minimum is determined over 
a single oscillation cycle (the use of a single cycle here 
is explained below). For N — 1500, the histogram is 
sharply peaked around the cell center at 1.0 ± 0.07 /im 
(1 standard deviation). As expected, with decreasing N 
the width increases: 0.09 /im at iV = 800, 0.16 /xm at 
N = 400, and 0.27 at iV = 200 The width of 

the midcell localization is large, particularly at protein 
cotmts {N — 200 and N — 400) that are significantly be- 
low those seen naturally. Hence, using significantly fewer 
protein copies degrades the accuracy of midcell division. 

Nucleoids, when present, also affect the positioning 
of the FtsZ ring through the poorly understood phe- 
nomenon of "nucleoid occlusion" where FtsZ rings 
do not nucleate over nucleoids and are restricted to ei- 
ther near the midcell or at the cell poles. Segregated nu- 
cleoids (at the 1/4 and 3/4 positions along the cell) will 
truncate the tails of the distributions shown in Fig. 3, 
further enhancing the accuracy of midcell division (in 
agreement with experiment [8ll26|). In normal cells with 
nucleoids, it is particularly important that the MinCDE 
system block polar FtsZ rings, since the nucleoids them- 
selves will inhibit FtsZ rings elsewhere away from mid- 
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FIG. 1: Ratio of average MinD density in right-hand 30% 
of the cell to that in left-hand 30%. Full line: stochastic dy- 
namics with random initial conditions at (a) A*' = 200 (b) 
N — 1500; dashed line: deterministic dynamics with equiva- 
lent parameters but with inhomogeneous initial conditions. 



cell. Assuming that FtsZ nucleation occurs at a single 
cycle MinD minimum (see below), then from Fig. 3 we 
see that N = 1500 is a high enough concentration to 
reduce the probability of polar division to considerably 
less than 0.01 per oscillation cycle. Given that about 50 
complete oscillation cycles normally occur between suc- 
cessive divisions, we see that attaining this level of ac- 
curacy is important. Significantly lower concentrations 
than N — 1500 will lead to an unacceptable probability 
of polar division, while higher concentrations will lead 
to only marginally increased accuracy, but at the cost 
of manufacturing many additional protein copies. From 
these simple arguments based on fluctuation effects, we 
see that E. coli may be using an optimal number of Min 
proteins, trading off midpoint precision against the cost 
of protein synthesis. There will be other constraints on 
the protein copy numbers (e.g., sufficient MinC to suc- 
cessfully inhibit off-center FtsZ ring formation), but fluc- 
tuations set useful boimds on the concentration levels. 

Comparison with experiment. Experimentally, the pre- 
cision of the MinCDE system can be probed in anucleate 
cells by measuring the position of the FtsZ ring. In these 
cells the only positional guide for division is the MinCDE 
system, which functions even in the absence of the nu- 
cleoids. Indeed in these anucleate cells the FtsZ ring po- 
sition is placed at midcell with a width of 0.12 /xm (scaled 
to a 2 /zm length) somewhat larger than the MinD 
distribution width 0.07 fim at TV = 1500. We would ex- 
pect that FtsZ ring nucleation would not precisely track 
the MinCD minimum, which could account for the dif- 
ference. In order to make this comparison we need to 
know how many cycles the FtsZ nucleation averages over 
to identify the location of the MinD minimum. Experi- 
mentally, this timescale is on the order of a minute, since, 
from Ref. 0, oscillation cycles of between 30 s and 120 s 
give normal division. As referred to earlier, this justifies a 
rough comparison of the width of the single cycle (~ 90 s) 
MinD density minimum distribution at iV = 1500 with 
that of the FtsZ distribution in anucleate cells. 
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FIG. 2: (a,b) MinD and (c,d) MinE relative density profiles in 
stochastic model with (a,c) = 200 (b,d) = 1500. Thick 
lines represent averages over (a,c) 160 and (b,d) 110 consecu- 
tive cycles respectively (maximum average density normalized 
to unity). Symbols (+, o, A, □) represent individual data sets 
that are averaged only over a single oscillation cycle. 
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FIG. 3: Histogram for position of MinD density minimum 
in the stochastic model with (a) N = 200, (b) N = 400, (c) 
N = 800, (d) N = 1500. The position of each minimum is 

computed by averaging over a single oscillation cycle. Data 
is for 160, 135, 120 and 110 consecutive cycles respectively. 



In conclusion, wc have studied fluctuations in the E. 
coli MinCDE system. In some regimes, we have found 
that fluctuations are essential for pattern generation. We 
have also found evidence that the MinCDE concentra- 
tions may be optimal for reliable midcell division. Based 
on the MinCDE system, we see that 0(10"*) copies of pat- 
tern forming proteins are required in bacteria to obtain 
positional information accurate to within a few percent. 
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